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Abstract. We present a series of hundreds of collisionless simulations of galaxy 
group mergers. These simulations are designed to test whether the properties of el- 
liptical galaxies - including the key fundamental plane scaling relation, morphology 
and kinematics - can be simultaneously reproduced by dry multiple mergers in galaxy 
groups. Preliminary results indicate that galaxy group mergers can produce elliptical 
remnants lying on a tilted fundamental plane, even without a central dissipational com- 
ponent from a starburst. This suggests that multiple mergers in groups are an alternate 
avenue for the formation of elliptical galaxies which could well dominate for luminous 
ellipticals. 



1. Introduction 

We aim to test the hypothesis that ellipticals form through hierarchical mergers of spi- 
ral galaxies in groups. If so, is the observed 'tilt' of the fundamental plane from the 
virial relation a natural consequence of hierarchical merging of spirals in ACDM cos- 
mology? Does the scaling relation of spirals - the Tully-Fisher relation - imprint itself 
in e lliptical properties as w ell? If so, this would contrast with evidence presented by 
e.g. [Robertson et al.l (1200 6:) . who suggest that elliptical properties (including the tilt of 



the fundamental plane) arise from varying gas fractions in binary mergers of spirals. 

Our method is to create a large statistical sample of galaxy merger simulations. 
These simulations naturally include mechanisms investigated independently in previ- 
ous works, including minor merging, hierarchical merging and group mergers. Similar 
simulations have been performed by lAceves & Velazquez! d2005l) . The chief improve- 
ment in this work is to use equilibrium, bar-stab le galaxy models based on M31 data 
dWidrow & Dubinskill200l IWidrow 6710 120081 and to fairly compare the results of 



hundreds of simulations to observations of local ellipticals (e.g. SDSS and Atlas3D). 



2. Simulations 

The initial conditions are designed to be quasi-cosmological - not an unbiased sam- 
ple but a roughly even sampling of the parameter space likely to produce groups with 
central ellipticals. Briefly, we create simulations using galaxy models with Sersic 
n s = 1 'pseudobulges' and n s = 4 classical/de Vaucouleurs bulges. Each model has 
2 sets of simulation with 7 target luminosity bins from 1/8 to 8 L* (in factors of 2), 
with 8 groups in each mass bins for a total of 2x7x8=108 simulations. Galaxies are 
placed within a sphere of r = r2oo of the group at z=l (or twice r2oo at z=l) with the 
most massive galaxy in the center. Each group has between N ga i s > N m i„ = 3 and 
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Figure 1. Images of a high-resolution re-simulation of a group of 10 galaxies 
totaling L*, starting from initial conditions and advancing 25,000 timesteps ( 5 Gyr) 
in each frame. 



Ng a i s < N max - 10 • sqrt(L/L*). Within each group luminosity bin, the number of 
galaxies in each simulation varies linearly from the minimum (3) to the maximum, so 
that L* groups have between 3 and 10 galaxies and the largest groups have 28 galax- 
ies. Galaxy luminosities are randomly sampled from an observed spiral luminosity 
(Schechter) function. Luminosities are drawn from a restricted range of the luminosity 
function with a width equal to (N ga i s + 2)/ 10 dex and such that the integral under the 
curve is equal to the target group luminosity. This setup excludes groups with a dom- 
inant bright galaxy and much smaller satellites (like the Milky Way and Magellanic 
clouds) which are unlikely to form ellipticals. A control sample of equal mass mergers 
is also created in each mass bin, one group using 3 galaxies and the other N max . Satel- 
lites are given preferentially inward and radial orbits, with the group as a whole being 

subvirial (to ensure collapse) and no satellite having |v| > v escape . 

Groups are simulated for 10 Gyr with the parallel tree code Partree (IDubinski 

19961) . using 50,000 timesteps and a softening length (spatial resolution) of 100 pc. 



Figure Q] shows a typical evolution for one such group. Simulations are analyzed after 
5, 7.5 and 10 Gyr. For brevity, only results for the 10 Gyr evolution (assuming groups 
formed at z=2) are shown. We note that groups with fewer galaxies typically merge 
in a few Gyr, while groups with more galaxies continue (slowly) accreting lower-mass 
satellites up to and beyond 10 Gyr. We process simulations with our own imaging 
pipeline, which is intended to create SDSS-like images of our groups, using sky back- 
grounds and S/N typical for SDSS imaging runs of galaxies at z=0.025. We also create 
kinematical maps at the same spatial resolution but with no noise in the velocity space. 
Each galaxy is imaged in ten equally spaced projections to boost our sample size by a 
factor of ten. 



3. Results 

Our main result is a measurement of the fundamental plane relation: log(R e ff) = 
a log(cr) + b/j + c. The tilt for our various samples is tabulated in table Q] Both pseudob- 
ulge and classical bulge mergers create a tilted fundamental plane. Figure |2] illustrates 
the small scatter and spatial extent of the plane. This challenges the interpretation of 
Hopkins et al.l (120081) that dissipation is necessary for the tilt, even if it is sufficient. 
Interestingly, equal mass mergers appear create a tilt. This could indicate that the tilt is 
a generic feature of multiple merging rather than being dependent on a particular spiral 
scaling relation; however, we caution that the equal-mass merger sample is quite small. 
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Table 1 . Best fit fundamental plane coefficients of central g alaxies in the simula- 
tions compared to SDSS observations (Hyde & Bernardi 2009). 



Sample 


N 


a 


b 


observed rms 


intrinsic rms 


Virial 


n/a 


2 


0.4 


n/a 


n/a 


Classical bulge 


101 


1.65 


0.296 


0.040 


n/a 


Exponential bulge 


85 


1.73 


0.327 


0.052 


n/a 


Classical bulge, equal mass galaxies 


33 


1.74 


0.291 


0.046 


n/a 


Exponential bulge, equal mass galaxies 


29 


1.72 


0.295 


0.054 


n/a 


SDSS 


50,000 


1.43 


0.315 


0.066 


0.058 


SDSS, stellar mass (excluding M*/L) 


50,000 


1.63 


0.336 


0.065 


0.049 



Fundamental plane Fundamental plane 




Figure 2. Face-on and edge-on projections of the fundamental plane of classical 
bulge mergers. 



The tilt appears similar to that from SDSS observations dHyde & Bernardill2009h 
excluding the M*/L contribution, which is nil by construction in our simulations. 

All galaxies in each image are fit with a single Sersic profile, although only central 
ellipticals are included in the final catalog. Classical bulges have higher central densi- 
ties than exponential bulges, so we would expect the inner bulge-dominated profile of 
classical bulge remnants to be steeper with larger n s . Equivalently, exponential bulges 
should be unable to form steep inne r profiles with large n s due to the conservation of 
phase space density dC arlber gl 1 1 98 6h . 

We confirm that classical bulge mergers have larger n s and a distribution more 
consistent with observations than exponential bulge mergers (Fig. [3]>. Classical bulge 
mergers appear to have slightly larger n s than the observations and exponential bulge 
mergers considerably lower. This suggests that ellipticals - if they are formed by dry 
mergers - are formed from spirals with a range of bulge profiles, consistent with mea- 
surements of spiral bulge profiles. 

Only classical group mergers produce a correlation between luminosity and n s . 
This relation is a result of more luminous ellipticals having been produced on average 
by m ore mergers. T he dependence of merger rate on halo mass is a prediction of ACDM 
(e.g. iHopkins et al.l ( 20101) ). However, exponential bulges are simply not concentrated 
enough to create merger remnants with n s > 4, even with repeated merging. Thus, 
luminous ellipticals are unlikely to be the product of exponential bulge mergers. 
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Figure 3. Sersic indices of central ellipticals. Classical bulge mergers (right) have 
larger n s for the same initial conditions and show a strong dependence of n s on galaxy 
luminosity, as with observed elliptica l s but u nlike exponential bulge mergers (left). 
Observed data from lNair & Abraham! (12010I) . 



We measure the traditional v/cr measure as the luminosity-w eighted average within 
R e ff, as used in IFU observations like Atlas3D (Cappellari et al.ll2011 ). As in Atlas3D, 
most of our ellipticals are slow rotators However, some remnants are formed with v/cr 
as large as 0.35, showing that it is possible to form fast rotators from dry mergers alone 
- a contrast with previous studies of dry binary mergers, which only formed very slow 
rotators. 

There does not appear to be any strong correlation between rotational support and 
number of mergers, galaxy luminosity or other parameters. We suspect the rotation is 
due to net angular momentum in the group. Whatever the case, it is possible for multiple 
mergers to produce fast rotators. However, fast-rotating SOs have been observed with 
v/cr > 0.5. If these SOs are merger products, they cannot have been formed from dry 
mergers alone. 
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